Exact solutions of the Saturable Discrete Nonlinear Schrodinger Equation 
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Exact solutions to a nonlinear Schrodinger lattice with a saturable nonlinearity are reported. For 
finite lattices we find two different standing-wave-like solutions, and for an infinite lattice we find 
a localized soliton-like solution. The existence requirements and stability of these solutions are 
discussed, and we find that our solutions are linearly stable in most cases. We also show that the 
effective Peierls-Nabarro barrier potential is nonzero thereby indicating that this discrete model is 
quite likely nonintegrable. 

PACS numbers: 61.25.Hq, 64.60.Cn, 64.75.+g 

The discrete nonlinear Schrodinger (DNLS) equation occurs ubiquitously [lj throughout modern science. Most 
notable is the role it plays in understanding the propagation of electromagnetic waves in glass fibers and other optical 
waveguides 2]. More recently it has been applied to describe Bose-Einstein condensates in optical lattices Here 
we are concerned with the DNLS equation with a saturable nonlinearity 



itp„ + (lpn+1 + 4>n-l - 2l/l n ) + 



1 + A*IVv, 



r^n = 0, 



(1) 



which is an established model for optical pulse propagation in various doped fibers fn Eq. (Q, ip n is a complex 
valued "wave function" at site n, while v and /i are real parameters. This equation represents a Hamiltonian system 
with: 
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so that Eq. is given by ivb n = J^-. The dynamics of Eq. conserve, in addition to the Hamiltonian TC, the 
power V 



N 



(3) 



X 



fn the above equations N is the number of lattice sites in the system. We note that a transformation \Jvvb n — > ip n 
will replace v by f and fi by ^ in the above equations. Note also that Eq. Q is invariant under the transformation 
ip n — > exp(i5)ip n where S represents an abitrary phase. 

For given system parameters v and /i it can be shown, using recently derived 5] local and cyclic identities for Jacobi 
elliptic functions 0, that Eq. Q has two (Case I and Case II) different temporally and spatially periodic solutions. 
Both solutions possess the temporal frequency 



w = 2 1 



v 
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Using standard notation 6] for the Jacobi elliptic functions of modulus m the solutions can be expressed as 
Case I: 



(4) 



, = 1 sn(/3,m) + c]p m) exp { _ i[ujt + d]) (5) 

y^icn(/3,m) 



where the modulus m must be chosen such that 



2m cn 2 (/3,m) 2K(m) 

v dn(/3,m) ' 1 N p ' 1 ' 
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FIG. 1: Illustration of the exact solutions of two types, v = 1, /i = 0.3, to = —1.33, and c = t — 5 = 0. 7V P = 5 (squares), 
N p — 10 (circles), and N p = 15 (triangles). Lines are guides to the eye. 



and c and S are arbitrary constants. We only need to consider c between and | (half the lattice spacing). Here 
K(m) denotes the complete elliptic integral of first kind [6(. While obtaining this solution, use has been made of the 
local identity 

it u , , \ i / m cn 2 (a,m) dn(a, m) . . , . 

dn (x, mjldmx + a, to) + an{x — a, to) = ttz -|dn(:r + a, to) + an{x — a, to) + 2 — — rdn(x, ml , (7) 

sn^(a,m) sn 2 (a,m) 

derived recently Q. In fact, given Eq. Q and this local identity [and similar ones for sn(x, m) and cn(x, to)], it was 
straightforward to obtain the two solutions presented here and the third solution follows simply by taking the limit 
to — > 1 of these two solutions as shown below. 
Case II: 



u /i dn(p, to) 



where modulus to now is determined such that 

2ii dn 2 Q3,TO) 4g(m) , Q , 

v cn(p, to) Aip 

While obtaining this solution, use has been made of the local identity [j| 

1/ sr / \ / m dn 2 (a, to). . . . cn(a,TO) 

men (x, mjlcnuc + a, to) + cn(x — a, m)J = — -|dn(x + a, to) + dn(x — a, m) + 2 — — rcn(x, to) . (10) 

sn^(a,m) sn z {a, to) 

Note that the two solutions, Eqs. JSJl and ©, are translationally invariant. 

The two solutions ipn 11 are illustrated in Fig. \I\for t = S = c = 0. In both cases the integer N p denotes the spatial 
period of the solutions. Both the solutions i\r n and ip^ 1 reduce to the same localized solution in the limit N p — > oo 
(m->l): 

Case III: 



^ = ±^p^l— e -^\ (N p ^oc), (11) 



*/7J cosh( [n + c]/3) " 

where /? is now given by 



sech/3=^. (12) 
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Again the frequency u> is given by Eq. J2J. This solution is noteworthy in that it is very similar in form to the 
celebrated exact soliton solutions of both the continuum cubic nonlinear Schrodinger equation and the (integrable) 
Ablowitz-Ladik lattice 

There are, as expressed by Eqs. J5J), JHJl, and l|12l) . stringent conditions on the parameters /i and v for which these 
exact solutions exist. In the cases I and II these limitations are illustrated in Fig. [2J which shows that the solution 
ipn only exists for parameter values below the lower curve (circles). Similarly, the solution tjj^ 1 for periods N p > 4 
only exists below the upper curve (squares). As can be easily seen from Eq. @ the ipn 1 solution does not exist 
for N p — 4. However, it does exist for N p = 3, but only for parameter ratios fi/is < 0. As a result of the periodic 
boundary conditions both solutions become meaningless for N p < 3. The solution tp^ 11 exists for all parameter values 
v > 2fi > 0. 




FIG. 2: Illustration of parameter values fi, v, and N p for which the exact solutions are allowed. Case I: 2fi/v between and 
cos 2 -j^- and N p > 3. Case II: 2^jl/v between and 1/ cos 2 j^- and N p > 3 except for N p = 4. 



For the ipn 1 solution, expressions for both the power Eq. and the Hamiltonian Eq. J2J) can be obtained by 
using exact (Poisson) summation rules 

V IU = l smh rg _ 2K(m)E(m) + 2K 2 {m)dn 2 (2K(m)c 1 m)l , (13) 
fx p l L J 



U UI = --sinh(/3) + ^] - Smh n }^ \(3 - 2K(m)E(m) + 2K 2 (m)dn(2K(m)c 1 m)] + -^2(3. (14) 
(i \ 2/i J /i p 2 - [i A 

Here the modulus m must be determined such that 

= *#^ sech^^, (15) 
A (mi) v 

where m\ — 1 — m is the complementary modulus and E(m) denotes the complete elliptic integral of the second kind. 
For the cases I and II analogous expressions can be obtained and they are given in the Appendix. 

In a discrete lattice there is an energy cost associated with moving a localized mode (such as a soliton or a breather) 
by a half lattice constant. This is called the Peierls-Nabarro (PN) barrier [Tol [ll| . Having obtained the expression 
for H 111 analytically in a closed form, we can now calculate the energy difference between the solutions when c = 
and c = 1/2, i.e. when the peak of the solution is centered on a lattice site and when it is centered half-way between 
two adjacent sites, respectively. We find that 

AE = n IH {c = 0) - H IH {c = 1/2) = sinh 2 (/3) sinh 2 (/3/2)A- 2 (m) < , (16) 

MP 



that is, the energy is lowest when the peak of the solution is centered at the sites. Thus, there is a finite energy 
barrier (i.e. the height of the effective PN barrier potential) between these two stationary states due to discreteness. 
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If the folklore of nonzero PN barrier being indicative of non-integrability of the discrete nonlinear system is correct, 
this suggests that quite likely our discrete model is non-integrable unlike the Ablowitz-Ladik model |g. 

In order to study the linear stability of the exact solutions ipi^ (j is I, II, or III) we introduce the following expansion 

i>n(t)=Tpi + tyn(t)e- iUt , (17) 

applied in a frame rotating with frequency u> of the solution. Substituting into Eq. Q and retaining only terms linear 
in the perturbation we get 

**. + <* +1 + - 2*,) + ( a + * + = » m 

Continuing by splitting the perturbation Stp n into real parts Su n and imaginary parts Sv n (Sip n = Su n + iSv n ) and 
introducing the two real vectors 

SU = {5u n } and SV = {Sv n } (19) 
and the two real matrices A = {A nm } and B = {B nm } by defining 



On,m-l + W — 2 H —77^ <W, (20) 




B n m = S n m+1 + <$n,m-l + W — 2 + " 7— ^nm, (21) 

V (i + mI^i 2 )/ 

where to ± 1 in the Kronecker 8 means: mil toog? iV. Then Eq. I|18l) can be written compactly as 

-SV + ASU = 0, and SU + BSV = 0, (22) 

where an overdot denotes time derivative. Combining these first order differential equations we get: 

SV + ABSV = 0, and SU + BASU = 0. (23) 

The two matrices A and B are symmetric and have real elements. However, since they do not commute AB 
and BA = (AB) T are not symmetric. AB and BA have the same eigenvalues, but different eigenvectors. The 
eigenvectors for each of the two matrices need not be orthogonal. 

The eigenvalue spectrum {7} of the matrices AB and BA determines the stability of the exact solutions. If it 
contains negative eigenvalues the solution is unstable. The eigenvalue spectrum always contains two eigenvalues which 
are zero. These eigenvalues correspond to the translational invariance (c) and to the invariance of the solution ip 3 n to 
a constant phase factor e~~ lS (i.e. translation in time), respectively In Fig. [3] we show the eigenvalue spectrum {7} 
for the cases I and II for several periodicities N p . It is important to note that in this figure we have N — N p . It 
turns out that the spectrum {7} is independent of c. The figure demonstrates that for N — N p , only the solution 
becomes unstable and this occurs only for N„ = 3. For all other values of N p both solutions are linearly stable. This 
also indicates that the localized solution V'n ^ s linearly stable; and we have checked that this indeed is the case in 
the entire existence interval. 

The solutions ip^, and ip^ 1 exist for all lattices N — JN p where J is a positive integer. However, we find ipn to be 
stable only for J = 1, while ip^ 1 is stable for all J. 

Finally, it is worth pointing out that Eq. also has an exact constant amplitude solution 

ip n (t) =ip exp[-i(ujt-qn + S)], (24) 

where S is a constant and u> satisfies the nonlinear dispersion relation 

c = 4sin 2 ( g /2)- "'^ |2 (25) 

where the wavenumber q = 2irp/N p in order to comply with the periodic boundary condition, and p is an intger. 

In conclusion, we have presented two spatially periodic and one spatially localized exact solutions of the DNLS 
equation with a saturable nonlinearity. We found these solutions to be linearly stable in most cases. We also 
calculated the Peierls-Nabarro barrier for the localized solution. These results are relevant for wave propagation in 
optical waveguides and doped fibers 0, , Bose-Einstein condensates as well as for many other nonlinear physical 
applications. Note that a related continuum version of Eq. (1), which arises in the context of the Fokker-Planck 
equation for a single mode laser, has been considered in Ref. [IJ- It would be important to search for ways of 
modifying the nonlinearity so that the PN barrier becomes zero-a possible route to an integrable model. 



This work was supported in part by the U.S. Department of Energy. 
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FIG. 3: Illustration of the stability of the exact solutions. Shown is the eigenvalue spectrum {7} for the matrix product AB, 
v = 1. Case I (left panel) and N p — 3 (triangles), N p = 4 (squares), N p = 5 (stars), and N p = 10 (circles). Case II (right 
panel) N p = 3 (triangles), N p — 5 (stars), and N p = 10 (circles). 



APPENDIX A 



In this appendix we give explicit expressions for H and P for the two spatially periodic solutions. While the 
importance of the energy expression is obvious, we would like to emphasize that the expressions for P could be 
used as a numerical diagnostic, for instance in keeping track of a conserved quantity in a simulation involving these 
solutions. 

Inserting the solution given by Eq. J5J into Eq. @ we get for the energy 

nI = 2 sn 2 (/?,m) / _ N ^ m) _ cs(/3; m)z(J3 m)) + f 1 _v\j- dn 2 ([n + c]/3, m)) 

+4 X>( 1 + 4^4^ + «1A m)) , (Al) 
p ~i V cn 2 (/3,m) ) 

where Z{(3 1 to) is the Jacobi zeta function and cs(/3, m) = cn(/3, m)/sn(/3, m). Also, use has been made of the identity 

dn(y, m)dn(y + a, m)—dn(a, to)— cs(a, m)Z(a, m)+cs(a, m)[Z(y + a, m)-Z(y, m)] and the fact that + 

1 + c), m) — Z(j3(n + c), m)\ = 0. From Eq. © we get for the power 

1 sn (p. to) \ 2 . l 

^ = 2^4^dn 2 ([n + c]/?,m). (A2) 

Similarly, inserting the solution given by Eq. JSJ into Eq. we get for the energy 



n 



2 sn 2 (/3,m) 
Mdn 2 (/3, m) 



N p (mcn(/3,m) - ds(0,rn)Z(J3,rn)) + ( 1 - — ) cn 2 ([n + c]/3, to) J 

V /V „ = i / 



A,. 



2 sn 2 (/3,m) 
A*dn 2 (/3, to) 



sn 2 (/3, to) 
dn 2 (/3, to) 

jv. 



cn 2 ([n + c]/3, to) 



AT p [TOcn(/3,TO) - ds(/3, to)Z(/3, to)] + f 1 - [ - (1 - m)N p + ^ dn 2 ([n + c]/3, to)] 

,1=1 

- +5> 

' 7 71=1 



cn 2 (/3, to) 
dn 2 (p\ to) , 



sn 2 (/3, to) 



1 + ~" 2 y n '"' dn 2 {[n + c]/3, to) 
cn 2 (p, to) 



(A3) 
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where again Z(f3,m) is the Jacobi zeta function and ds(/3, to) = dn(/3, m)/sn(/3, m). Also, use has been made of the 
identity |5j m cn(y, m)cn((/ + o, m)=m cn(a, m) - ds(a, to)Z(o, m)+ds(a, to) [Z(y + a, m)-Z(y, m)]. From Eq. (J2Jl we get 
for the power 

„, r 1 sn 2 (/5,m) v-A 9 , r ,„ , 1 sn 2 (/?,m) / , . v-^ 9/r . 1 . 

V U = - 2 7 E cn 2 ([n + c]/3, m) = - - »' -JV p (l - m + E ^ w + c ^ w ) ' ( A4 ) 
A 4 dn (/3, to) ^ M dn (/3, to) I ^ / 

In order to get the sums over the same expressions for and V 11 as for and V 1 we have used the basic relations 
cn 2 (x, to) + sn 2 (ir, m) = 1 and dn 2 (ir, m) + TOsn 2 (ir, to) = 1. In the continuum limit (small (3, large N p ) the sums may 
be replaced by integrals. First 

^ d „ 2fl „ +c]Am) ,^ = «^)|W =Wp |W, (A5) 

' n n Kim) Kim) 

n— 1 v y v ' 

where Q = 2 in Case I and Q = 4 in Case II. The other sum 

^, / sn 2 (/3,m) , 2/r .„ A / 7t8 2 (/3,to) \ 

z — ' 1 cn z (p,TO) / \2vl — TOA (m)cn^(p, to)/ 



n=l 



where 6(/3, to) is the Jacobi theta function. For to — > 1, Eqs. (Al) and (A3) can be used to determine the asymptotic 
interaction between two nonlinear solutions given by Eq. (11). 
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